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ABSTRACT 

Using a Gibbs sampling algorithm for joint CMB estimation and component separation, we compute 
the large-scale CMB and foreground posteriors of the 3-yr WMAP temperature data. Our parametric 
data model includes the cosmological CMB signal and instrumental noise, a single power law fore- 
ground component with free amplitude and spectral index for each pixel, a thermal dust template with 
a single free overall amplitude, and free monopoles and dipoles at each frequency. This simple model 
yields a surprisingly good fit to the data over the full frequency range from 23 to 94 GHz. We obtain 
a new estimate of the CMB sky signal and power spectrum, and a new foreground model, including a 
measurement of the effective spectral index over the high-latitude sky. A particularly significant result 
is the detection of a common spurious offset in all frequency bands of ~ —13 /iK, as well as a dipolc 
in the V-band data. Correcting for these is essential when determining the effective spectral index of 
the foregrounds. We find that our new foreground model is in good agreement with template-based 
model presented by the WMAP team, but not with their MEM reconstruction. We believe the latter 
may be at least partially compromised by the residual offsets and dipoles in the data. Fortunately, the 
CMB power spectrum is not significantly affected by these issues, as our new spectrum is in excellent 
agreement with that published by the WMAP team. The corresponding cosmological parameters are 
also virtually unchanged. 

Subject headings: cosmic microwave background — cosmology: observations — methods: numerical 



1. INTRODUCTION 

A major challenge in CMB research is component sepa- 
ration, which can be summarized in two questions. First, 
how can we separate reliably the valuable cosmological 
signal from confusing foreground emission? Second, how 
can we propagate accurately the errors induced by this 
process through to the final analysis products, such as 
the CMB power spectrum and cosmological parameters? 

During the last few years, a new analysis framework 
capable of addressing these issues in a statistically con- 
sistent approach has been developed. This framework 
is Bayesian in nature, and depends critically on the 
Gibbs sampling algorithm as its main comp utational en- 
gine. The pioneering ideas were d escribed bv lJewell et alJ 
(|2004h and lWandelt et all (|2004f ). and later implemented 
for modern C MB data set s for te mpera ture and po- 
larizat ion by lEriksen et alJ ([2004a) and lLarson et alJ 
(|2007l ). respective ly. Applications to the 1-yr and 3- 
yr WMAP data Bennett et al.l l2003al ; iHinshaw et al 



120071: [Page et alj|2007l) w ere described by lO'Dwver et al 
( 2004L and lEriksen et all (j2006t l2007al ). These paper; 



papers 

mainly focused on the cosmological CMB signal, and 
adopted the foreground corrected data provided by the 
WMAP team. 
Recently this algorithm was extended to include inter- 
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nal com ponent separation capabilities by lEriksen et alJ 
(|2007blh Using very general parameterizations of the 
foreground components, this method produces the full 
joint and exact foreground-CMB posterior, and there- 
fore allows us both to estimate each component sepa- 
rately through marginalized statistics and to propagate 
the foreground uncertainties through to the final CMB 
products. The implementation of this algorithm used in 
this paper is called "Commande r" , and is a d i rect de - 
scendant of the code presented by lEriksen et al.l (|2004al) . 

In this Letter, we apply t he method to the 3-yr WMAP 
temperature observations (|Hinshaw et al.l l2007f ). This 
data set, with five frequency bands, allows only very 
limited foreground models; however, the analysis pro- 
vides a powerful demonstration of the capabilities of the 
method. For a comprehensive analysis of a controlled 
simulation wi t h ident ical properties to this data set, see 
lEriksen et aD (|2007bl ). 

2. DATA 

We consider the 3-yr WMAP temperature data, pro- 
vided on Lambda 7 in the form of sky maps from ten 
"differencing assemblies" covering the frequency range 
between 23 and 94 GHz. Since our current implemen- 
tation of the Gibbs foreground sampler ca n only handle 
sky m aps with identical beam response (|Eriksen et al.l 
2007b), we downgrade each of these maps to a common 
resolution of 3° FWHM and repixelize at a HEALPix 8 
resolution of iV S ide = 64, corresponding to a pixel size 
of 55'. These ten maps are then co-added by frequency 
into five single frequency band maps at 23, 33, 41, 64 
and 94 GHz (K, Ka, Q, V and W-bands, respectively). 

The power from the instrumental noise is less than 1% 

7 http://lambda.gsfc.nasa.gov 

8 http://hcalpix.jpl.nasa.gov 
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of the CMB signal at £ = 50 in the V- and W- bands, and 
less than 2% at I = 100 (jEriksen et al.ll2007bf ). To regu- 
larize the noise covariance matrix at high spatial frequen- 
cies, we added 2 /iK per 3° pixel of uniform white noise. 
This noise is insignificant at low multipoles, but domi- 
nates the signal near the spherical harmonic truncation 
limit of £ max = 150. We then have five frequency maps at 
a common resolution of 3° FWHM, with signal-to-noise 
ratio of unity at £ ~ 120, and strongly bandwidth limited 



at 



150. 



We choose to include such high Vs in the analysis for 
two reasons. First, our main goal is an accurate approx- 
imation of the CMB likelihood at I < 50. In order to 
ensure that the degradation process (i.e., smoothing and 
noise addition) does not significantly affect these multi- 
poles, it is necessary to go well beyond I ~ 80-100. Sec- 
ond, significant information on the spatial distribution of 
foregrounds is obtained by going to higher resolution. 

The cost of this treatme nt of the noise is high x 2 val- 
ues in the ecliptic plane (jEriksen et alJ l2~007bl ). where 
the instrumental noise is higher because of WMAP's 
scanning strategy. However, since these high x 2 val- 
ues are caused by unmodelled smoothed, random, white 
noise, they do not indicate a short-coming of the signal 
model, but only a slight under-estimation of the statis- 
tical errors on small angular scales. This has been con- 
firmed by otherwise identical analyses at both 4° and 
6° FWHM. We present the 3° FHWM case here as a 
compromise between angular resolution and accuracy of 
the noise model. This issue will further suppressed with 
additional years of WMAP observations and, eventually, 
high-sensitivity Planck maps. 

We impose the base WMAP Kp2 sky cut 
([Bennett et al.1 l2003bh on the data, but not the 
point source cuts. The base mask is downgraded from 
its native iV S idc = 512 resolution to N S ia c = 64 by 
excluding all low-resolution pixels for which any one of 
its sub-pixels is excluded by the high-resolution mask. A 
total of 42 081 pixels are included, or 85.6% of the sky. 

The frequency bandpass of each map is modelled as a 
top-hat function, and implemented in terms of effective 
frequency as a funct ion of spectral index as described by 
lEriksen et al.l (|2006h . The freque ncy specifications of t he 
WMAP radiometers are given bv lJarosik et al.l j2003). 

3. MODEL AND METHODS 

We adopt the following simple parametric model T v {p) 
for the observed signal (measured in thermodynamic 
temperatures) at frequency v and pixel p, 



T v (p) = s{p) + m° v + ™>t [e< ' n(p)] 



+b 



t(p)a(v) 



.dust 



1.7 



f(jp)a{v) [ — j 



v \ 



0{p) 



(1) 



The first term is the cosmic CMB signal, character- 
ized by a frequency-independent spectrum and a co- 
variance matrix in spherical harmonic space given by 
the power spectrum, {a* tm atm') = Ce5a'5 m m>, where 
s (p) — J2e m a ^ m (p)- The secon< i an d third terms 
denote a free monopole and three dipole amplitudes at 
each frequency. We use the standard Cartesian basis 



vectors projected on the sky, {l,x, y, z}, as basis func- 
tions for these four modes. The fourth term represents 
a template-based dust model scaled by a spectral index 
of P = 1.7, in which t(p) is the dust template (FDS) of 
iFinkbeiner et al.l (|1999l ). evaluated at i^ ust = 94 GHz, 
and a(t) is the conversion factor between antenna and 
thermodynamic temperatures. The last term is a sin- 
gle foreground component with a free amplitude f(p) 
and spectral index (3(p) at each pixel. The reference fre- 
quency for this component is — 23 GHz. 

The free parameters in this model are: 1) spherical har- 
monic coefficients ag m of the CMB amplitude s; 2) CMB 
power spectrum coefficients Cg; 3) monopole and dipole 
amplitudes at each band; 4) the amplitude of the dust 
template; and 5) amplitudes and spectral indices of the 
pixel foreground component. 

The WMAP data do not have sufficient power to con- 
strain this simple completely by themselves, as there is a 
very strong degeneracy between the foreground compo- 
nent amplitudes at each pixel and the free monopole an d 
dipole coefficients at each band (jEriksen et alJ [2007bf ). 
For this reason, we introduce two priors in addition to 
the Jeff ereys' ignorance prior discussed bv lEriksen et al.l 
(20073). First, we impose a Gaussian prior on the spec- 
tral indices of (3 3 ± 0.3 : A d irect fit of the 408 

MHz template (jHaslam et al.l 119811 ) to the WMAP K- 
band data for Kp2 sky coverage imp l ies an index of -3 
(jDavies et al.ll2006l ), and lDavies et al.l (j!996f ) determined 
a typical range for high latitude spectral indices between 
408 MHz and 1420 MHz of -2.8 to -3.2. Note that this 
prior has a noticeable effect only at high Galactic lati- 
tudes, where the absolute foreground amplitude is low. 
At low galactic latitudes, the data dominate the prior 
by up to a factor of ~ 50, and any potential bias in the 
near-plane free-free regions is negligible. 

Second, we impose an implicit spectral index orthogo- 
nality prior o n the monopo l e and d ipole coefficients, as 
described bv lEriksen et alJ (|2007bl ). projecting out the 
frequency component of these coefficients that matches 
the free spectral index map, thus effectively determining 
the zero-level of the foreground amplitude map. We also 
tried an alternative approach, first estimating the Q, V, 
and W-band monopole and dipole coefficients separately 
given a crude estimate of the spectral index map, and 
then estimating all other parameters given these coeffi- 
cients. Results were very similar. Thus, the two priors 
adopted in this analysis have a very weak effect on all 
main results. 

Having defined our model and priors, we map out the 
joint posterior distrib ution using th e foregr ound Gibbs 
sampler described by lEriksen et al.l (j2007bh . We refer 
the interested reader to that paper for full details of the 
algorithm, and for a comprehensive analysis of a realistic 
simulation corresponding to the same data and model 
used in this Letter. 

Finally, we estimate a new set of cosmological parame- 
ters within the standard A CDM model. For this analysis, 
we follow the approach of lEriksen et al.l (|2007af ). and re- 
place the low-£ part of the WMAP like lihood with a new 
Blackwell-Rao Gibbs-based estimator (jChu et al.l [2005). 
No ancillary data sets beyond the 3-yr WMAP temper- 
ature and polarizatio n data are included in the analysis. 
The CosmoMC code (jLewis fc Bridlell2002D is used as the 
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Fig. 1. — Marginal posterior mean maps in Galactic coordinates. 
Rows from top to bottom show the CMB reconstruction, the fore- 
ground amplitude, and the foreground spectral index, respectively. 



TABLE 1 

MONOPOLE AND DIPOLE POSTERIOR STATISTICS 





Monopole 


Dipole X 


Dipole Y 


Dipole Z 


Band 




(A*K) 


QxK) 


(UK) 


K-band 


-11.8 ±0.5 


1.7 ±1.2 


-2.2 ±0.8 


2.2 ±0.1 


Ka-band 


-16.6 ±0.5 


0.9 ± 1.2 


0.9 ±0.8 


-1.3 ±0.1 


Q-band 


-12.8 ±0.5 


1.9 ± 1.2 


-0.9 ±0.8 


0.4 ±0.1 


V-band 


-11.1 ± 0.5 


1.6 ± 1.2 


-3.9 ±0.8 


4.0 ±0.1 


W-band 


-12.6 ± 0.5 


1.7 ±1.2 


-0.9 ±0.8 


1.0 ±0.1 



Note. - Means and standard deviations of the 
marginal monopole and dipole posteriors. 

main MCMC engine. 

4. RESULTS 

Figure Q] shows the marginal posterior mean maps for 
the CMB sky signal, the foreground amplitude and the 
foreground spectral index. Table [T] gives the correspond- 
ing results for the monopole and dipole coefficients for 
each frequency band. The FDS dust template ampli- 
tude relative to 94 GHz and an assumed spectral index 
of = 1.7 is b = 0.917 ± 0.003. The CMB power spec- 
trum is discussed separately below. 

Figure[2]shows the average x 2 computed for each Gibbs 
sample. A \ 2 value exceeding x 2 = 15 corresponds to 




Fig. 2. — Mean x 2 map computed over posterior samples in 
Galactic coordinates. A value of \ 2 = 15 is high at the 99% sig- 
nificance level. 



rejection of the model in that pixel at 99% statistical 
significance. Two features are clearly visible in this plot. 
First, the ecliptic plane, or rather, WMAP's scanning 
strategy, is clearly visible, and this is mainly due to the 
unmodelled smoothed noise component at high fs, as 
discussed in Section [2] 

Second, there are clearly visible structures near the 
Galactic plane, and in particular around regions with 
known high free- free emission (e.g., the Orion and Ophi- 
ucus regions). This is likely due to the fact that a single 
power-law is not a good approximation to the sum of 
many foregrounds with comparable amplitudes. Thus, 
we have clear evidence of modelling errors in this so- 
lution, and we therefore strongly emphasize that the 
quoted error bars presented in this paper include for- 
mal statistical errors only, and not systematic, model- 
dependent irncert^amties. These issues are discussed in 
depth by lEriksen et al.l (|2007bl ). who find similar be- 
haviour for a controlled, simulated data set. 

Although x 2 is high, the CMB solution obtained is ev- 
idently good. First of all, in the ecliptic poles, where the 
WMAP instrumental noise is suppressed by the scanning 
strategy, the x 2 distribution is essentially perfect. This 
implies that the signal model as such is adequate at high 
latitudes. Second, the CMB map is virtually without 
signatures of residual foregrounds. (Note that the signal 
inside the Galactic plane is partially reconstructed using 
high- latitude information and the assumption of isotropy. 
The signal on scales smaller than the mask size is washed 
out because it is not possible to predict these from higher 
latitudes.) Finally, the foreground amplitude and spec- 
tral index maps correlate very well with known templates 
of synchrotron and free-free emission. For instance, the 
spectral index near the Gum nebula and near the Vela 
regions are close to (3 = —2.1, as expected for free- free 
emission. 

The single most surprising aspects of the solution are 
the monopole and, possibly, the V-band dipole coeffi- 
cients, listed in Table Q] Most notably, there is a strong 
detection of a roughly 13/iK offset common to all fre- 
quency bands. Formally speaking, these offset values 
are only optimal within the current model; however, this 
type of signal is not degenerate with any other compo- 
nent in the model. Further, we have attempted to fit 
several other models assuming no offsets at one or more 
bands. These all result in strong, visible residuals in 
the CMB map, and considerably higher x 2 values over- 
all. Finally, very similar results have been obtained by 
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Fig. 3. — Difference maps of the total "Commander" W-band 
foreground model with the WMAP MEM model (top panel) and 
with the WMAP template fit model (bottom panel). 

other researchers 9 through other methods, although, to 
our knowledge, these results have not yet been published 
in the literature. We therefore believe that the monopole 
and dipole coefficients presented here are more optimal 
even in an absolute sense than those obtained by WMAP 
based on a cosecant fi t to a plane parallel galaxy model 
([Bennett et al.ll2003bl ). 

Figure [3] shows the difference maps between the 
Commander W-band foreground model and the MEM 
and te mplate-based foreground models of lHinshaw et alJ 
(2007). Clearly, our foreground model agrees surpris- 
ingly well with the simple template fits, but not with the 
MEM solution. One possible expl anation for this i s that 
although the MEM approach of iHinshaw et al.1 (|2007l ) 
does attempt to estimate spectral indices for each pixel, 
it does not include monopole or dipole components in its 
model. Therefore, the MEM solution is plausibly com- 
promised by the non-zero offset detected here, at least 
in part. In addition, the MEM method could be biased 
by the initial subtraction of the (foreground contami- 
nated) ILC estimate of the CMB anisotropy from the 
frequency maps, and the use of the 408 MHz data as 
a prior. Conversely, Commander could be compromised 
to some extent by the use of power law spectral indices 
for the combined low frequency foreground component. 
Nevertheless, the difference is surprising given that the 
W-band foreground is expected to be mostly comprised 
of thermal dust emission. 

The marginal maximum posterior CMB power spec- 
trum is shown in Figure |4j together with the maximum- 
likelihood/ pseudo-Cf hybrid spectrum computed by the 
WMAP team. Perhaps the most notable difference is 
in the I = 21 multipole, which looks anomalous in the 
WMAP spectrum, as noted by other authors. 

The cosmological parameters corresponding to the 
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Fig. 4. — The CMB temperature power spectra obtained in 
this paper (red) and that by t he WMAP team (bla ck). The best- 
fit ACDM model spectrum of[Spcrg;cl ct al. (2007]) is shown as a 
dashed line. 

Commander spectrum for a standard six-parameter 
ACDM model are tt b h 2 = 0.0222±0.0007, ft m = 0.243± 
0.036, log(10 10 A s ) = 3.027 ± 0.068, h = 0.730 ± 0.032 
and r = 0.089 ± 0.030. Corresponding values for the un- 
modified WMAP likelihood are n h h 2 = 0.0221 ± 0.0007, 
O m = 0.242 ± 0.035, log(10 10 A s ) = 3.030 ± 0.068, 
h = 0.730 ± 0.032 and r = 0.091 ± 0.030. These val- 
ues refer to marginal means and standard deviations. 

Clearly, the agreement between the two sets of results 
is excellent, and this provides a strong confirmation of the 
WMAP results: At the level of precision of the WMAP 
experiment, details in the foreground model used for fore- 
ground correction appear to have only a minor impact on 
the CMB temperature power spectrum. 

5. CONCLUSIONS 

We have presented the first exact Bayesian joint 
foreground-CMB analysis of the 3-yr WMAP data. We 
have established a new estimate of both the CMB sky 
signal and the power spectrum, a detailed foreground 
model consisting of a foreground amplitude and spectral 
index map and a dust template amplitude, and also pro- 
vided new estimates of the residual monopole and dipole 
coefficients in the WMAP data. 

The detection of significant non-zero offsets in the 
WMAP data is the new result of the greatest immedi- 
ate importance for the CMB community. These new 
monopole and dipole estimates could have a significant 
impact on several previously published results, espe- 
cially those concerning the foreground composition in 
the WMAP data. For example, our foreground model 
is in excellen t agreement wi t h the simple template fits 
presented by Hins haw et al.l (|2007l ). but not with their 
MEM reconstruction. 

Taking a longer perspective, the most important as- 
pect of this analysis is a demonstration of feasibility of 
exact and joint foreground-CMB analysis. This will be 
essential for Planck, whose high sensitivity and angular 
resolution demand more accurate foreground separation 
than WMAP. Considering the flexibility, power, and ac- 
curacy of the method employed in this paper, together 
with its unique capabilities for propagating uncertainties 
accurately all the way from the postulated foreground 
model to cosmological parameters, we believe that this 
should be the baseline analysis strategy for Planck on 
large angular scales, say I < 200. 
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All results presented in this paper and the basic Gibbs 



samples are available at http://www.astro.uio.no/~hke 
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